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; Abstract 

O ■ Low-temperature (< 60K) phase changes of free clusters containing 50, 137, and 229 

■ CH4 molecules have been observed in isoenergetic Molecular Dynamics computations. Bulk 
solid methane exhibits structural phase transformation at 20.4 K. However, clusters of 50 

Q^' molecules already melt at ~ 25 K, which makes the observation of solid-solid transformations 

rather challenging. 

^ ! 
> 

^ 1 Introduction 

o ; 

^ ■ Small free clusters of molecules with very high symmetry (octahedral) have been found to resemble 
^ bulk of the same substance even in the case of small number of ingredients (30 to 80) pp. The 
O . number of particles (molecules or/and atoms) that mimics bulk properties like phase transitions 
^ ! is strongly dependent on the range of the potential. The effect of decreasing the range of the 
\ potential is to destabilize strained structures [2^. In the case of long-range potentials (Coulomb 
[T]) one detects structures, determined by the co-operative (collective) interaction - the system 
'-^ • is expected to transform easier than in the case of short-range interactions. On the other side, 
. the surface of free clusters plays much more important role for long-range potentials. Various 
• ^ parameters of the potential, for instance electric charge value, govern the cluster response to 
^ \ external changes. An example is the melting temperature increase for higher charge values 
d In the present work we study clusters consisting of less symmetrical molecules (tetrahedral) in 
order to check how the lower symmetry of the ingredients influences cluster's behavior at different 
temperatures. The interacting potential contains a long-range term (Coulomb), which makes 
it necessary to study the size influence - clusters of 50, 137, 229 methane molecules have been 
simulated with the help of the isoenergetic Molecular Dynamics method described in the next 
section. Simulations of methane clusters have been performed more than two decades ago 
They showed that the lower limit of the classic approach is ~10K. In those days the computers 
were less fast and the researchers tended to go for cheap rather than reliable potentials. Nowadays 
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more sophisticated potentials are available and better results could be obtained This 
justifies re-visiting the problem of how the cluster size infiuences the structural phase changes. 
Another intriguing question is how a free cluster, e.g. a cluster at zero external pressure, behaves 
when the temperature is rather close to the ultimate low limit of the classic mechanics calculations. 

Solid (bulk) methane exhibits several crystalline phases based mainly on the difference in the 
orientation of molecules in the unit cell. Infrared investigation of solid methane performed by Bini 
and co-workers jH] let him propose a phase diagram, which has been additionally elucidate by 
Nakahata et al. with a help of optical and X-ray diffraction studies [Zj . The liquid methane freezes 
at 90.6 K and p = 0.1 MPa [HI into a fee orientational disordered phase I JU], characterized with 
rotational diffusion of the molecules. A transition to the orientational ordered phase II occurs at 
20.4 K JU] at zero pressure. The phase II is the James-Keenan cubic structure, in which one out 
of four molecules (Oh) is undergoing a weakly hindered rotation, while the other three molecules 
{D2d) have an anti-ferro order [TT] . 

The state of orientational disorder takes place for the systems of spherically- symmetric 
molecules {CH^ with a symmetry lower than the site in the crystal structure 0. At cooling such 
systems could undergo phase transitions to a partially-ordered or fully ordered state. Hindered 
rotations occur because of weakly angular-dependent intermolecular forces and a large rotational 
energy as it is in methane [TT] . 

Is such a behavior inherent for small systems as well? This is the question we study in the 
current work. 



2 The model 

Due to their spherical symmetry methane molecules resemble in many aspects the noble gases (Ar, 
Kr, Xe). The intra-molecular frequencies (~ 4 x 10^^, ~ 9 x 10^^ Hz) are one order of magnitude 
higher than the frequency of the intermolecular vibrations. That is why, we consider the molecules 
as rigid bodies (C-H distance b=1.094 A). Suitable intermolecular potentials for regarding molec- 
ular space orientations could be the Lenard- Jones potential jHj or a more sophisticated 3-body 
RMK potential ^21- In both approaches complicated procedures had been followed to obtain pa- 
rameters' values. Needless to say, these parameters are not entirely satisfactory. In our study, the 
total potential is a sum of pair atom-atom potentials Upw{i,j): 
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where a , /3 denote either a carbon or a hydrogen atom; is the distance between the i-th and 
the j-th atom; gj „ is either qc = —0.572 or = 0.143 in elementary charge units. 
The parameters have been obtained in a way to fit experimental data. For instance, the charge 
and the bond-length are determined from the octopole moment value, aap in A, eap in meV sue 
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Figure 1: Lindemann index for clusters of 50 
and 229 molecules; a 6cc-starting configuration. 



Figure 2: Lindemann index for clusters of 50 
and 229 molecules; a /cc-starting configuration. 



0: 

acc = 3, 35 , auH = 2, 81 , och = 3.08 ; 
ecc = 4.257 , enn = 0.715 , ecn = 1, 7447. 

Two starting configurations with randomly oriented molecules have been designed: 1) a bcc lat- 
tice, which resembles the most spherically-symmetric shape of a cluster and 2) fee, which is a 
possible lattice for a bulk methane. We execute the velocity Verlet algorithm E] for numerical 
integration of classic equations of motion in a special micro-canonical ensemble - N, E, p=0 (free 
clusters) at temperatures above lOK to avoid quantum effects. The integration step of dt = Ifs 
guarantees an accuracy greater than 0.003% of the energy conservation for runs up to ~ 0.5ns; 
the records are taken every 100 /s. 



3 Results 

Liquid- and solid-like phases can be distinguished with the help of a specially constructed formula, 
which is an extension of the Lindemann index 6iin ^1] : 



N{N - 1) 



N 
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where Vij is the distance between molecular centers of mass. In bulk, the Lindemann index is 
computed on the basis of the mean deviation of the i-th atom from its ideal lattice position and 
6iin < 0.1 indicates a solid state. In a free cluster the surface plays a destabilizing role and the 
cluster is solid if < 0.08 [T3] . 

The Lindemann index, plotted in the Figs.(P - Ej), shows three distinguishable phases : two 
solid-like phases below 30irand a liquid-like phase above 30i<'for all three sizes and both starting 
configurations. 

The caloric curves (total energy per particle as a function of temperature) are given in the FigEl 
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Figure 3: Lindemann index for heating-cooling- 
heating of 137 molecule cluster; a 6cc-starting 
configuration. 



Figure 4: No hysteresis is observed on 'cooling- 
heating' for the two starting configurations; 
only the case of a fecc-starting configuration is 
shown. 



and Figini Although the caloric curves are too smooth (a change of the slope is hardly seen, no 
more than 6%), we observe a region of frustration in the temperature interval (18,25) K. In our 
computations the average temperature of the system is computed from its average kinetic energy: 
< Ej^in >= 3/2fc^T, with ks - the Boltzmann constant. Analyzing the trajectories at different 
total energies, we see that the system temperature jumps unevenly in the above interval. This is 
an indication for changes much more clearly seen with the help of the Lindemann index. There is 
no hysteresis on 'heating-cooling' of the system - FiglU 

The different phase structures are distinguished on the basis of their radial distribution func- 
tions - Eq. 121 the radial distribution of molecular centers of mass shows the structure of the 
cluster at a specific temperature, while the atom-atom radial distribution reveals orientational 
order (disorder) of the molecules. 

V ^ 

^w = ]^(EE'^(^-^'^)) (2) 

i=l j^i 

Fig El shows the existence of a solid- like structure below 30 and the absence of that structure 
above 30 K. The volume of the cluster enlarges at heating as well. Fig |H1 shows typical distribution 
of the atoms for the oriented (solid line in the figure) and disoriented (dashed line in the figure) 
solid-like phases. The first peak of the distribution correspond to the H-H neighbors, the second 
peak corresponds to the C-H neighbors, and the third one is for the C-C neighbors. The same 
maximum of the H-H and C-H peaks indicates a unique mutual orientation of the molecules, e.g 
the cluster is orientationally ordered. The phase transition at ~ 20K causes disorientation of the 
molecules as it is seen in the Fig|Hl- the dotted curve. 

Finally, we present the Density of States (DOS) function obtained by a normal mode analysis 
of the quenched system ^H!- The FigO shows evidently the existence of a rovibrational coupling 
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Figure 5: Caloric curves at heating of clusters Figure 6: Caloric curves for the same cluster 
containing 50, 137, and 229 molecules, a bcc- sizes with a /cc-starting configuration, 
starting configuration. 
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Figure 7: Radial distribution of the molecular 
centers of mass for a 50 Ci74 cluster; a /cc- 
starting configuration. 



Figure 8: Atom-atom radial distribution for 50 
CH4 cluster, /cc-start. 
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Figure 9: The DOS for a 229 CH^ cluster, a 
/cc-starting configuration at 11 K - the total 
spectrum is indicated with 'all' in the figure 
inlet; 'rot' is the rotational and 'vib' is the vi- 
brational spectrum. 



Figure 10: Vibrational spectrum for 229 Cif4 
cluster; fee - start. 



1^]. The FiglTUl shows broadening of the vibrational spectrum at heating. 

4 Conclusion 

Implementing Molecular Dynamics method to a microcanonical state of free methane clusters, 
we show that they exhibit transformations between an orientationally-ordered (at ~ lOK) and 
disordered phase (above 20K). The lack of hysteresis in the cooling- heating regime indicates a 
continuous transition. However, such a statement can be confirmed either by studying infinite 
systems (i.e. periodic boundary conditions) or by applying the finite-size scaling theory [T7] . 
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